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ABSTRACT 


Terminal voltage measurements with various cathodes and 
anodes in a high power, quasi-steady MPD discharge show that 
the magnitude of the current at the onset of voltage fluctua- 
tions is an increasing function of cathode area and a weaker 
decreasing function of anode area. Tests with a fluted cath- 
ode indicate that the fluctuations originate in the plasma 
adjacent to the cathode rather than at the cathode surface. 

Measurements of radiative output from an optical cavity 
aligned to examine the current-carrying portion of a two- 
dimensional, 56 kA MRD discharge reveal no lasing in that 
region, consistent with calculations of electron excitation 
and resonance radiation trapping. A voltage-swept double 
probe technique allows single-shot determination of electron 
temperature and electron number density in the recombining 
MPD exhaust flow. 

Current distributions within the cavity of MPD hollow 
cathodes for various static prefills with no injected mass 
flow are found to replicate distributions measured with in- 
jected flows, yielding the important result that the dynamics 
of the injected flow does not play an important role in the 
hollow cathode emission process. 
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I . INTRODUCTION 


i 

During the preceding six-month interval, technical papers 
summarizing our progress in all three areas of investigation, 
MPD discharges, plasmadynamic lasers and hollow cathodes, were 
presented at AIBA meetings. In the MPD portion of the program, 
a summary of the work on current-limiting phenomena was pre- 
sented at the AIAA 11th Electric Propulsion Conference, March 
19-21, 1975 at New Orleans, Louisiana. One graduate student 
is now working in this area of perfoinnance limitations in MPD 
accelerators. His recent work, included in this report, shows 
additional encouraging results for MPD devices using argon as 
a propellant. 

In the plasmadynamic laser studies, our program of two 
complementary approaches continues. The principal results of 
the first of these, the study of the collisional-radiative 
recoinbination mechanism in a three-dimensional MPD exhaust, 
were reviewed in a paper presented at the AIAA 8th Fluid and 
Plasma Dynamics Conference, June 16-18, 1975 at Hartford, 
Connecticut. More recent work in this program has led to the 
development of a single shot electron temperature measurement 
technique to aid in determining the electron cooling rates in 
the recombining flow. The goal of the second plasmadynamic 
laser project is the direct demonstration of lasing in a 
unique two-dimensional discharge configuration. In this 
report, the results of a spectroscopic study of the current 
carrying portion of the discharge with a resonant optical 
cavity are detailed. 

Like the MPD discharge work, the principal results of the 
hollow cathode program were presented at the AIAA Electric 
Propulsion Meeting at New Orleans. The approach outlined in 
that paper, namely the examination of the active zone length 
at low currents and mass flows, has been extended to even 
lower currents. The recent addition of a second graduate 
student in this program allows a broader examination of cur- 
rent conduction mechanisms in hollow cathodes. 
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IX. QUaSI-STEAEY MPD DISCHaRGE 
Cathode Studies (Rudolph) 

The terminal voltage of an MPD accelerator has been found 

to exhibit highly oscillatory behavior for discharge currents 

above a certain value. The onset of these oscillations corre- 

lates with the parameter, J /m, where J is the onset current 

at which the oscillations arise and m is the propellant mass 

flow.^ 1/166 although earlier investigations have shown that 
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this parameter is a strong function of the cathode geometry, 
the previous report showed preliminary data indicating that 
J ^/m also depends weakly on the anode geometry. ^ , Recent 
results support this anode area dependence. The results show 
that for a 5.0-cm-long, stainless steel cathode, the onset 
current increases by less than a factor of two as the anode 
area is decreased by a factor of ten. For larger cathodes the 
effect becomes even weaker. 

Earlier floating probe measurements on several cathodes 

revealed that the voltage fluctuations originate in the cathode 

166 

region of the discharge. Unfortunately, these studies were 

unable d distinguish between a cathode surface effect (related 
to the electron emission process) and a volume effect in the 
plasma surrounding the cathode. To differentiate between these 
two possibilities, the performance of a fluted cathode was 
investigated, and the results, presented herein, indicate a 
cathode plasma effect is responsible for the oscillations. 

A. Discharge Configurations 

Figure 1 shows the two anode geometries and examples of 

various cathode geometries used in this study. The small anode 

2 

with a surface area of 250 cm is shown in Fig. la. it con- 
sists of a 1-cm-thick aluminum plate, with an inner surface 
facing the 5.1-cm-deep discharge chamber, a lip region with a 
minimum diameter of 10,2 cm, and an outer face with a maximum 
diameter of 18.6 cm. The outer barrel of the anode, which is 


26.3 cm 


GEOMETRIES 

FIGURE I 

AP25- 5119 
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approximately 34 cm long, is electrically insulated in this 

2 

geometry. The larger anode has an area of about 2500 cm and 
is shown in Figs, lb, c, and d. It is identical to the small 
anode, with the exception that the outer barrel is no longer 
insulated. This larger anode has been used in all previous 
investigations of the voltage fluctuations. 

The cathodes were all stainless steel cylinders of various 
lengths, each with a hemispherical tip. The radius was identi- 
cal for all cathodes in order to maintain a constant geometric 

61 

terra in the electromagnetic thrust equation, while the length 
was changed to vary the cathode area. 

B. Onset Current 

Figure 2 shows a graph of the onset current versus cathode 

area for the large and small anodes. The lower curve, measured 

166 

with the large anode, has been discussed in previous reportsT 
The upper curve represents the results of the study using the 
small anode. For all of these data, the mass flow of argon 
propellant was 6 g/sec. 

The results using the small anode are qualitatively simi- 
lar to those using the large anode except that the entire curve 
has been displaced upwards. This relatively weak dependence of 
critical current on anode area is reasonable insofar as changes 
in the anode can be expected to affect the current and poten- 
tial profiles throughout the discharge, including the region 
surrounding the cathode. Nevertheless, the inverse depend- 
ence of onset current on anode area presents a logical problem 
if taken to the extreme of zero anode area. Without further 
data, it is anticipated that additional current limiting 
phenomena would become manifest for veiY small anode areas. 

It should be noted that the abscissa in Fig. 2 equivalent- 
ly represents the cathode length since all data were obtained 
with cathodes of the same radius, i.e. these data do not allow 
discriminating between a cathode area effect and a cathode 
length effect. 
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Another iiiteresting implication of the data is that the 
largest value of the critical current found experimentally 
(30 3^A) gives a value of J /ih of 150 kA *sec/g. Using this 
value in the electromagnetic thrust equation yields a specific 
impulse of 2500 sec, which if achieved in practice provides a 

significant improvement over the 1250 sec value measured hy 

168 

Boyle and the 850 sec limit originally proposed by 
Malliaris.^”^ 

Most of the previous data obtained during the study of 
the onset current dependence on cathode area was done using 
tungsten cathodes. 0?he switch to stainless steel was made due 
to the relatively low cost and ease of machining. In addition, 
because the work functions of the two materials are comparable, 
it was felt that there would be no difference in the discharge 
characteristics. To confim this expectation experimentally, 
a tungsten cathode was used with the small anode and the 
result is shown by the data point labelled W in Fig. 2 . The 
critical current for the tungsten cathode falls on the curve 
generated by the stainless steel cathodes, thereby confirming 
the interchangeability of tungsten and stainless steel. 

C. Origin of Potential Fluctuations 

While the previous floating probe measurements indicated 
that the voltage fluctuations originate in the cathode region, 
they did not allow differentiation between a cathode emission 
process and a plasmadynamic effect centered about the cathode. 
An experiment was subsequently designed to aid in this dis- 
crimination since a cathode emission mechanism would result in 
a surface area dependence while an effect originating in a 
plasma layer separated from the cathode surface would be 
independent of the fine details of the surface. 

Figure 3 shows a stainless steel cathode with 12 flutes 
machined into its surface. The flutes were not carried to the 
upstream edge of the cathode in order to leave the injection 
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geometry (which depends on the base of the cathode as a channel 

wall) unchanged. The 12 flutes are each 0.3 cm deep with a 

spacing between flutes of 0.2 cm. The surface area of this 

fluted cathode is approximately 48 cm , as compared to a sur- 

2 

face area of 30 cm for the previously studied cylindrical 
cathode of the same outer radius and length. 

Estimates of the Debye length in the cathode region of the 
discharge give a value of 10 cm. This shows that the emis- 
sion process takes place essentially at the cathode surface, 

even within the 0 . 3 cm deep flutes . The electron— ion mean free 

—2 

path is estimated at 10 cm, confirming the penetration of the 
plasma into the regions between the flutes. Any azimuthal var- 
iations in electron number density introduced by the fluted 
cathode should be attenuated over a radial distance of approxi- 
mately 10 ^ cm since the random thermal velocity of the elec- 
trons is roughly 100 times their radial drift velocity. Thus 
the plasma conditions surrounding this cathode should be virtu- 
ally unchanged from that of the standard cylindrical cathode at 
the same current. This implies that if the oscillatory process 
originates in the plasma surrounding the cathode, the critical 
current should be approximately the same for both the fluted 
'athode and the cylindrical cathode, whereas if the fluctua- 
tions originate at the cathode surface, the critical current 
for the fluted cathode may be up to 50% greater than that of 
the cylindrical cathode. 

The results of the fluted cathode test, shown in Pig. 4, 
indicate that the critical current for the fluted cathode is 
the same as that for the cylindrical cathode. The arguments 
above then lead to the conclusion that the oscillatory process 
is based on an as yet unidentified plasma phenomena in the 
region just off the cathode surface and is independent of the 
cathode emission process. Having thus localized the source of 
the fluctuations, it follows that all of the onset current 
data to date, including that of the fluted cathode, can be 
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reduced to a single graph (for a fixed anode area) by using 
cathode length rather than cathode surface area as the 
independent variable. 

D. Summary 

Since the previous report, several experiments have been 
conducted to improve our understanding of the physical nature 
and origin of the terminal voltage fluctuations which arise 
under certain operating conditions in the MPD arcjet. An un- 
expected result is that the onset current depends on the anode 
geometry as well as the cathode geometry, in the inverse direc- 
tion, although the anode dependence is much weaker. Experi- 
ments with a fluted cathode have shown that the fluctuations 
originate in the plasma surrounding the cathode rather than 
at the cathode surface itself. Further studies are being pur- 
sued to delineate more sharply the region from which the 
oscillations first arise. 
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III. PiASMaDyNaMic laser studies 

A. MPD Discharge in a Laser Cavity (Dutt) 

Previous semi-annual reports'^ ' have followed the pro- 
gress of a set of experiments intended to determine directly 
whether lasing can be sustained in the discharge and exhaust 
regions of an MPD accelerator. The discharge apparatus is a 
two-dimensional version of the conventional MPD accelerator 
with parallel plates as electrodes. A direct evaluation of 
lasing is possible by placing the plasma region under consider- 
ation within a resonant optical cavity. The use of a cavity 
avoids the experimental difficulties sometimes associated with 
the direct measurement . of the small signal gain. Furthermore, 
the two-dimensional geometry, if .sustaining a uniform discharge 
can reduce to a minimum the effects of variation of plasma 
properties and emission coefficient along the optical axis of 
the cavity. 

This report briefly describes the most recent configura- 
tion of the apparatus, and some of the operational character- 
istics of the discharge. A search for lasing was conducted in 
some parts of the discharge and the results are also summarized 

1. Resoiiant Cavity and Discharge Apparatus 

The resonant optical cavity remained as described in the 

1 66 

previous semi-annual report. It consists of 2 concave 

mirrors of reflectivity 99.7% and 97% in the wavelength range 
0.4 to 0-53 \x. Their separation is slightly larger than their 
radius of curvature, 147 cm. The cavity can be aligned 
externally using a He-Ee laser. 

The discharge apparatus is a 45-cm-wide parallel plate 
configuration consisting of two 10-cm' long anodes on either 
side of a 5-cm-long cathode, all made from aluminum. In its 
present form, the lateral edges of the electrodes are recessed 
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into 40-cm-high by 56-cm-long insulating side plates coated, 
with boron nitride as shown schematically in Pig. 5. In an 
earlier version of this device/ these sidewalls were very 
small (10 cm X 5 cm) and extended only as far as the elec- 
trodes , 

2 . Experimental Results 

a. Discharge Characteristics 

The discharge characteristics of this configuration were 
established earlier in the development of this program. The 
main concern was the achievement of a laterally uniform current 
distribution within the discharge. This is a prerecjuisite for 
lateral uniformity of the emission coefficient across the dis- 
charge. It has been demonstrated that for this electrode con- 
figuration, the discharge occupies the entire width of the 
electrodes only at currents above 40 kA. For these conditions, 
with the small insulating sidewalls, measurements using a mag- 
netic probe indicated that the current density was a maximum 
near the center of the electrodes and decreased towards the 
edges. In the present discharge apparatus with the larger 
insulating sidewalls, the current distribution is more uniform 
than before, as shown in Pig. 6. 

Another useful indicator of discharge behavior is the 
terminal voltage. Oscillograms of terminal voltage V were 
recorded for various total currents J and mass flows m. 

The voltage vjas found to be quasi-steady in most cases, a 
typical trace of which is shown in Pig. 7 for a mass flow of 
20 g/sec. 

The variation of the quasi-steady voltage with current is 
shown in Pig. 8 for a mass flow of 11.4 g/sec, and the varia- 
tion with mass flow is shown in Pig. 9 for a current of 56 kA. 
An extrapolation of the curve in Fig. 8 to J = 9 indicates 

that V{J = 0) = V = 50 V at 11.4 g/sec. V has been iden- 

o o 

tified as the sum of the cathode and anode fall voltages. 



FIGURE 5 

AP 25-5133 





PLEXIGLAS 


ALUMINUM 


SIDEWALL WITH 


VIEWING PORTS 


T\A/0 DIMENSIONAL DISCHARGE APPARATUS 






MAGNETIC FIELD DISTRIBUTION 

J = 56 k A , y = 2 cm , Z = 5. 1 cm 


15 


J-71440 



I, , 34KA/D1V 

UPPER ANODE CURRENT 


100 V/DIV 

TERMINAL VOLTAGE 


sec 



J-7l^^ 4b 

ilHIHMHI 

!SS!S I2.34 kA/DIV 


LOWER ANODE CURREr.'T 


TERMINAL VOLTAGE AND CURRENT 
J = 56 kA,m=ECg/sec 


FIGURE 7 

AP 25- P-566 



VOLTS 


16 



CURRENT, kA 

CURRENT-VOLTAGE CHARACTERiSTtC 


FIGURE 8 

AP20 - 5140 




FIGURE 9 

AP25-9I4I 


500 


400 


300 




o 

> 


200 


100 


0 




18 


163 

present even if the interelectrode space has zero resistance. 
Log-log plots of V - V^vs J are straight lines as shown in 
Fig. 10 for mass flows of 11.4 and 25.9 g/sec, and can there- 
fore he represented in the form 

(V - Vq) = k (1) 

where the exponent n is dependent on the mass flow. 

Magnetic field and terminal voltage measurements provide 
general information about the discharge characteristics. In 
order to gain further insight into the species distribution in 
succeeding levels of ionization, spectroscopic diagnostics 
were employed. 

b. Spectroscopic Study 

The insulating sidewalls of the discharge apparatus have 
openings with removable plugs which permit optical sighting of 
the plasma transverse to the discharge (along the x direction 
in Fig. 5) . Using these openings, spectra can be recorded at 
several locations within the current conducting region and in 
the exhaust plume. In addition, the resonant cavity can be 
positioned to enclose these same locations, and spectra can be 
recorded through its output mirror, which transmits 3% of the 
light. 

If lasing occurs, then the radiation intensity in the 
laser lines is much greater than in the non-lasing transitions. 
This would be immediately obvious in spectra photographed 
through the output mirror of the cavity. A spectrum photo- 
graphed at the same location but without the cavity would pro- 
duce a very different distribution of line radiances. There- 
fore the event of lasing can unambiguously be demonstrated by 
this comparison of emitted spectra. 

To aid in the selection of suitable transitions where 
lasing may be expected, it is useful to review briefly the 
processes that can lead to overpopulation of the upper level 
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of a laser transition. One such process is the excitation 

mechanism postulated for the operation of conventional argon 
A“-2 

ion lasers. In these devices, the upper level is populated 
by inelastic collisions between energetic electrons and argon 
atoms in their ground state. ^Dhis process may proceed by 
either a one-step excitation directly to the upper laser level 
or by a two-step excitation with the ion ground level serving 
as an intermediate state. Theoretical estimates and some ex- 
periments indicate that the one-step excitation is significant 
only for current densities ^ 500 A/cm , assuming an electron 
temperature sufficiently high to overcome the threshold energy 
for this process (35.5 eV) . At current densities in the 
neighborhood of the threshold current density for lasing, 
50-70 A/cm , the two-step process is dominant, but even this 
mechanism requires a mean electron energy in excess of about 
3 eV to sustain an adequate population inversion in conven- 

tional ion lasers Because the current density anticipated 

. . 2 
xn the two-dimensxonal device is on the order of 100 A/cm , 

the two-step excitation process may be capable of producing 
inversions in the current carrying portion of the MPD dis- 
charge. 

An alternate process for creating inversions is the 
recombination of a dense plasma of low electron temperature. 

In this case, the inversion is fed from higher lying levels 
of the same species or by recombination from the next higher 
state of ionization. Although lasing in argon as a result 
of pumping by this process has not yet been successfully 
demonstrated, estimates of the inversion ratio made in this 
laboratory indicate that the requisite low electron tempera- 
ture and the high electron density exist within the expanding 
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plume of the MPD exhaust flow. 

The configuration of the discharge apparatus permits 
spectroscopic diagnosis of both the discharge region, where 
elevated electron temperatures may create inversions by 





excitation from lower levels, and the exhaust fltjw where popu- 
lation inversion by recombination may predominate. During the 
period covered by this report, the possibility of lasing was 
investigated only within the discharge region. 

Spectra were photographed both with and without the reso- 
nant cavity. The locations where spectra were taken through 
the resonant cavity are designated by their y, z coordinates as 
(0, 0.95 cm); (0, 3.81); (0, 10.8); (1.80, 6.90); (2.92, 2.16) 
and (4.05, 9.13). Reference spectra, photographed without the 
cavity, were also recorded at some of these locations. The 
y and z coordinates of the optical axis are defined in Pig. 5 
where the origin of the coordinate system is the center of the 
front edge of the cathode-, 

The study was confined to a current level of 56 kA where 
the discharge is known to distribute itself uniformly across 
the width of the electrodes. The other variable in the control 
of the experimenter, the injected mass flow, was varied from 
0.55 g/sec to 55 g/sec. 

Two spectra, typical of those taken without the cavity, 
are shown in Pig. 11 for mass flows of 20 and 54 g/sec. Most 
of the lines in these spectra, which cover the range from 
0.45 to 0.5 |a, correspond to 4s - 4p transitions in the argon 
ion. All. Some weak lines of AI in this spectral region do 
not appear in Pig. 11, but since these lines are not generally 
prominent, a definite statement of the relative abundance of 
AI and All cannot be made from these spectrograms - 

Spectra taken through the output mirror of the cavity 
exhibit (at a considerably reduced intensity) about the same 
intensity distribution between different lines and between the 
different mass flow rates as observed without the cavity. It 
must therefore be concluded that no lasing v;as observed for 
these discharge conditions at these locations. 
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3. Discussion 

One explanation for the absence of lasing in the dis- 
charge region may be the excitation rates, which are based in 
turn on cross sections for electron impact excitation of ions. 
In order to correlate the results of the present investiga- 
tions with conventional argon ion laser operations, the cur- 
rent density and the ratio of electric field to pressure, E/p, 
can be compared. It has been previously stated that the cur- 
rent densities calculated from magnetic field measurements 
are typically 50 to 100 A/cm in the MPD discharge, as com- 
pared to threshold current densities for lasing of 50 to 70 
A/cm in conventional ion laser operation. The B/p parameter 
has been found useful in correlating argon ion laser opera- 
tion and is therefore used here to compare the present work 
with known laser conditions. The pressure within the dis- 
charge region was estimated from known flow rates to be 0.3 
to 5 torr, a range which overlaps typical ion laser pressures 
of 0.05 to 0,5 torr. However, the electric fields estimated 
from terminal voltage measurements indicate that for this dis- 
charge configuration E/p is usually less than 0.2 V/cm*torr, 
which is significantly lower than the E/p values observed in 
ion lasers of about 25 V/cm*torr, 

The energy gained by an electron in an electric field is 
proportional both to E and to the mean free path between 
energy absorbing collisions. Since the mean free path is in- 
versely proportional to gas pressure, the electron energy is 
proportional to E/p, Therefore, the low E/p in the present 
experiments may imply an electron temperature too small to 
sustain a significant population inversion. 

The lower E/p in the MPD device results, in the first 
instance, from the low voltage in the MPD arc compared to the 
higher voltage of the glow discharges in conventional ion 
lasers. Figure 9 shows that it is possible to increase the 
terminal voltage in the MPD discharge by simply decreasing the 
mass flow. However, even the modest increases in E field that 
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might he realized by this procedure are counteracted by the 
measured trend for the current contours to bow further down- 
stream as mass flow decreases, thus decreasing the E field by 
increasing the current path length. 

A second possible cause of the absence of a population 
inversion is resonance trapping of the lower level of the 
transition. The lower level of the 0.4880 p All laser transi- 
tion normally depopulates by radiative decay to the ion ground 
state. This resonance transition has a high probability, 

A, = 2.64 X 10 /sec, when compared to that of the laser tran- 

7 

sition itself, — 6.59 x 10 /sec, a disparity which helps 
to maintain a population inversion even when excitation to the 
lower level is comparable to excitation to the upp^r level . 
However, if there are a large number of particles in the 
ground state, then the resonance radiation may become trapped, 
leading to a decreased effective probability of decay for the 
lower level and a consequent reduction or elimination of the 
inversion. The detrimental effects of resonance trapping can 
be expected to increase with plasma dimensions since the 
greater the distance that the resonance radiation has to 
travel before leaving the plasma, the greater is its probabil- 
ity of reabsorption and trapping. 

Effects of resonance trapping in our discharge configura- 
tion can be expected to be severe compared to conventional ion 
lasers. The ion laser, usually maintaining the discharge in a 
thin tube, is of very small lateral extent ( < 1 cm typically) 
and there is consequently little resonance trapping in the 
active volume. The large dimensions (20 cm minimum) of the 
present discharge configuration precludes easy escape of 
resonance radiation before reabsorption. 

In conclusion, a combination of low electron temperature 
and excessive resonance trapping offers a possible explanation 
for the absence of lasing in the current conducting region of 
the MPD discharge. Similar spectroscopic studies are presently 
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in progress in the expanding region of the MPD exhaust, where 
the collisional-radiative recombination model appears to be 
more promising for producing an inversion, 

B. Electron Temperature Measurement (Campbell) 

The goal of the plasmadynamic laser program is the under- 
standing of the relaxation processes which occur in the exhauet 
pltime of a high power, quasi-steady MPD arc. In a sufficiently 
rapid expansion, i.e, one in which the electron temperature 
relaxes in a distance short compared to the recombination 
length of the plasma, these processes can lead to population 
inversions between bound electronic states of the accelerated 
species. Although theoretical studies for hydrogen have 
specified a range of electron density and temperature for which 
lasing is possible by this collisional-radiative recombination 
mechanism, similar calculations have not been performed for 
argon due to its complex internal structure. It is thus of 
interest to determine experimentally the electron temperature 
profile in the exhaust flow, from which the electron cooling 
rate can be deduced. 

In the past the electron temperature in the quasi-steady 
MPD flow has been measured using a double Langmuir probe tech- 
nique. The voltage-current characteristics for these probes 
were constructed point -by-point by applying a fixed voltage to 
the probe electrodes from a small capacitor bias circuit and 
then measuring the probe current during the discharge. 

The determination of the electron temperature at a single loca- 
tion thus required many tens of shots in order to guarantee 
probe reproducibility and to establish the requisite accuracy 
of the characteristic. 

For a complete mapping of the electron temperature dis- 
tribution such as is planned in the plasmadynamic laser pro- 
gram, a new technique has been developed in which the voltage 
bias on the double probe is swept at a preselected time and 
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rate during the discharge, Q?his method, discussed below, 
allows a complete characteristic to be generated v/ith a single 
firing of the arc, 

1. Probe and Circuit Description 

The double probe used in this study consists of two egual- 
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area tungsten wires, 7.6 x 10 cm in diameter, mounted in a 
quartz support. The wires have an aspect ratio of 100 to 
avoid end effects, and have a separation of 100 diameters to 
insure negligible interaction between the two electrodes. 

A special voltage biasing unit has been developed for 
these measurements capable of providing either a free-running 
sawtooth wavefom or just 1 to 4 segments of this waveform. 

The amplitude can be varied from -20 to +20 volts, the sweep 
time per segment can be as short as 50 psec (compared to a 
discharge time of % msec) , and in the partial segment mode, 
the time of voltage sweep initiation can be delayed with 
respect to an initial timing mark . 

A schematic illustration of the double probe circuitry is 
shown in Fig. 12. Probe current was monitored with a Tektronix 
type 6042 current probe, and probe bias was measured differen- 
tially using two Tektronix P-6013 (1000 x) voltage probes. 

The probe voltage sweep was triggered by a Hewlett-Packard 
pulse generator which was delayed several hundred microseconds 
with respect to discharge initiation to allow the development 
of steady plasma conditions at the probe location. For 
most of the data reported here, only one segment of the saw- 
tooth (a linear ramp) was used, i.e. the bias circuit main- 
tained a constant negative bias on the probes until triggered, 
at which time it swept to a bias of equal amplitude but oppo- 
site sign. 

2. Data Reduction 

Using the Debye length Aq , the smallest mean- free -path 
A , and the probe radius R, it is convenient to fom the 
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dimensionless parameters \/R and X/R whicli serve to define 
the regimes where an electric probe can operate. Using values 
of the electron temperature and number density expected in the 
exhaust plume, upper and lower bounds of these ratios are; 

8.7 X 10"“^ ;i^/R -i 8.7 X 10"^ and 0.25 < A /r -t 40. 

Prom the bounds of /R, it can be concluded that the 

probe will be operating in the thin sheath regime. However, 
because the Knudsen number, X /R# be of the order 1, the 
shadowing effect of the probe must be considered, i.e., the 
probe does not subtend a negligible solid angle at the point of 
the last charged particle collision before absorption by the 
probe. When severe, this problem is manifest by the lack of a 
linear relation between the log of the electron current and the 
probe bias in the electron retarding region of the characteris- 
tic, thus prohibiting the determination of electron temperature 
in the conventional manner. 

Talbot, et al, have shown that an accurate value for T 

e 

can still be obtained from a double probe when operating in a 

A— 7 

transition regime. From their results, the electron temper- 
ature can be found from the following relation 


(1 +‘ 3 ‘)| ( 2 ) 

where A^ and A^ are the areas of the two probes, jj^ and are 

the ion saturation currents, ^ V/^l is the inverse slope of 

the characteristic, and CT is a correction factor dependent on 

Xip /R, and the probe aspect ratio. This relation is 

essentially the collisionless result obtained by Johnson and 
A “*8 

Malter using the equivalent resistance method, modified by 
the correction factor 1 + cT , For the entire range of plasma 
parameters expected in the plume, CT is of the order of lO”^, 
and thus the classical theory can be used while incurring an 
error of only a few percent. 


T (^K) = 11,600 
0 




v=o 
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The prohe axis for all tests was aligned with the local 
velocity vector. The resulting V - J characteristic is iden- 
tical to that which would he obtained in a stationary plasma 
provided that 


i_ (k 

Ao u 


(3) 


where £ is the probe length, the ion mass, and u the 
plasma velocity A- 10 present case, V - & (10^) 

verifying that the characteristic measured in the flowing 
plasma is valid. 


Figure 13a shows a record of the probe current (upper 
trace) and the arc current (lower trace) . The probe bias for 
this case was fixed at a value sufficiently large to produce 
ion current saturation. It can be seen that following the 
passage of the initial plasma front, the probe current is 
steady for approximately 500 lisec. Therefore, it was decided 
to trigger the voltage sweep approximately 600 psec after dis- 
charge initiation. Figure 13b shows an expanded trace of the 
arc current, with the region where the probe characteristic 
was measured highlighted by increased luminosity. The probe 
current, taken for a fixed bias on a 20- times expanded time 
scale during the highlighted portion of Fig, 13b and displayed 
in Fig. 13c, shows that the ion saturation current is steady 
over the interval in which the voltage is swept. 

Figure 14 shows a typical record of the probe voltage 
and current over a 100 psec interval during which the bias 
was swept from -4 to +4 volts. Because the voltage increases 
linearly in time, the lower trace in Fig. 14 represents the 
cur rent-volt age characteristic of the double probe multiplied 
by a scale factor which is constant and depends on the sweep 
rate. Thus the electron temperature, can be obtained from 
Eqn. 2 (without the small correction) , using the saturated 
current densities from the oscillogram and calculating the 
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slope of the characteristic from, 

dV 
dt 

This procedure was used for all of the data reported in the 
next section. 

3. Results of Probe Measurements 

The experiments reported here were designed to verify that 
electron temperature measurements using the single shot, swept 
voltage technique are reliable. For all tests, the probe was 
mounted on the centerline, 32 cm downstream of the anode, and 
the arc current and mass flow were maintained at 4 ka and 
12 g/sec respectively. This operating condition and probe 
location were chosen because previous optical depth measure- 
ments indicated a population inversion among the states of the 
0.4880 \x argon ion transition for this case. 

One of the principal concerns in determining a temperature 

by a voltage sweep technique is that the bias change will be 

too rapid to obtain the correct value of the temperature. 

Calculations show that the highest frequency associated with a 

bias change of the order of kT^/e (8 x 10 sec ) is small 

compared to the ion plasma frequency evaluated at the electron 

8 —1 

temperature (7 x 10 sec ) , i*e. the probe sheath is easily 

A— 11 

capable of adjusting to the changing probe bias. Still, 

the ultimate test of the data lies in comparing the deduced 
temperatures over the maximum range of biasing rates. 

To test the effect of frequency, the voltage sweep time 

and amplitude were varied from 50 to 400 lasec and 5 to 20 volts 

4 5 

respectively, r rresponding to 2.5 x 10 xi dV/dt<8-5xlO V/sec. 
Over this entire range of sweep rates, the measured electron 
temperature was the same within 5% at a value of 2760 °K. 

As an aside, this electron temperature measurement may 
help to explain the population inversion that has been indicated 
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by other data. The low temperature not only rules out colli- 
sional excitation from lower lying states, but also is neces- 
sary to explain the large recombination rates which result in 
a flux of electrons from the continuum into the bound states 
of All to produce the inversion. 

Further evidence that the plasma frequency response is not 
imposing any limitation on the acquisition of single shot tem- 
peratures is provided by the probe current- volt age characteris- 
tics in Fig. 15. These data were obtained using a voltage 
sweep rate of 8.5 x 10^ V/sec with a probe whose aspect ratio 
was 200. In this figure, the solid line represents the charac- 
teristic determined by plotting the instantaneous values of the 
swept probe bias and current from a single oscillogram; the 
superimposed data points were obtained by setting a fixed bias 
on the probe and measuring the probe current during the 
500-|isec quasi-steady phase of the discharge. Each data point 
represents the mean of several shots using the same bias . The 
fact that the curves are virtually identical is conclusive 
proof that the voltage sweep technique is reliable. 

In the above experiments the probes were glow-cleaned 
after every firing of the arc to ensure that the electrode sur- 
faces were not contaminated, thus producing a spurious current- 
voltage characteristic. In order to determine the influence of 
material deposited on the probes, characteristics were obtained 
from successive firings of the accelerator without cleaning the 
probes. Figures 16a, b and c show the probe current after 1, 4 
and 7 discharges. Comparing Figs. 16a and b, the first effect 
observed is a reduction in the slope of the probe current at 
V = 0, which produces an incorrectly large {dV/dI)^_Q. 

Since the saturation currents aren't appreciably altered, use 
of this characteristic would result in an erroneously large 

value of electron temperature. This effect has been observed 
A— 12 

by others and is apparently caused by charges in the probe 
contact potential due to the deposition of thin layers of for- 
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eign material. After many arc firings, the characteristic is 
considerably distorted (Pig. 16c) . The hint of two saturation 
regions may be due to self-cleaning of the probe at higher 
voltages, exposing more of the electrode surface area. 

In summary, the voltage sweep technique offers a rapid 
and reliable method for determining the current-voltage charac- 
teristic of a double probe. With this diagnostic, it is now 
possible to obtain spatial maps of electron temperature over 
a wide range of arc operation. These mappings have been initi- 
ated with the goal of a better understanding of the relaxation 
processes occurring in the exhaust flow of the accelerator. 



37 


IV. HDIIiOW CATHODE STUDIES (Krishnan) 

in the preceding semi-annual report, a comprehensive pic- 

X66 

ture of the Princeton hollow cathode program was presented. 

In it were described experimental measurements of current and 

potential distributions in various large hollow cathodes for 

—3 

arc currents from 0.9 to 7 and argon mass flows from 10 
to 16 g/sec. These measurements showed that at the highest 
current the current distribution within the cavity is uninflu- 
enced by changes in the cathode configuration, whereas varying 
the current and mass flow for a fixed cathode configuration 
produces significant changes in the current conduction pattern. 
For a given current, maximum current penetration into the cav- 
ity occurs at mass flows of approximately 10*”^ g/sec, decreas- 
ing for both higher mass flows, up to 10 g/sec, and lower mass 
flows, down to 10 g/sec. For a fixed mass flow, the penetra- 
tion monotonically increases as the current decreases. Spec- 
troscopic measurements of AI and All line radiation, and spec- 
tral photographs taken through a 4880 A All filter confirmed 
these attachment patterns. 

In the most recent work, the current and potential distri- 
butions were measured at a still lower current of 0.25 kA for 
the range of argon mass flows from 10 to 0.4 g/sec. When 
these results are compared with earlier results obtained with 
the same mass flows but with higher currents of 0,9 to 7 kA, a 
consistent picture emerges of the effects of current and mass 
flow on the current distribution within a large hollow cathode. 

A. Experiments 

To determine first whether the maximum penetration of cur- 
rent into the cavity continues to increase as the current is 
decreased below 0,9 kA, the current distributions were measured 
for a current of 0.25 3cA and several argon mass flows. 



38 


The current distributions were obtained by traversing a 
0.3-cm-dia, magnetic field probe at a fixed radial position 
of 0,8 cm. Using standard passive integration techniques and 
the measured asimuthal symmetry of the cavity plasma, the 
probe directly yields the current enclosed in a circular cross 
section whose radius equals the probe radial position. Assum*- 
ing purely radial current flow between the O.S*’cro radial posi- 
tion and the cavity wall (E. = 0.95 cm), the local slope of the 
enclosed current profile yields the surface current density. 

£*igure 17 shows a typical enclosed current profile and 
its associated radial current density profile for HC XII, 
shown for comparison at the top of the figure. The arc oper- 
ating conditions for these data are a total current of 0.25 kA 
and a mass flow of 0.1 g/sec argon. As a check on the current 
density data, integration of the profile over the cavity sur- 
face gives a total current of 0.254 kA, well within the esti- 
mated 10% error bar for these data. 

Figure 18 shows the surface current density profile from 
Fig. 17, along with several other profiles taken at various 
mass flows from 5 x lO" to 0.4 g/sec. The current was fixed 
at 0.25 kA and the cathode was again HC XXI, From these dis- 
tributions of surface current density, three characteristic 
features emerge: 

1) As the argon mass flow is reduced from 0,4 g/sec to 
0.1 g/sec, the peak in the surface current density moves from 
0.3 cm to 1.9 cm upstream of the end of the cathode cavity, 

2) At the same time, the current attachment at the sur- 
face becomes more diffuse, leading to a drop in the peak 
current density. 

3) Further reduction in the mass flow from 10~^ g/sec 
to 5 X 10 g/sec causes the peak in. the current density dis- 
tribution to move downstream towards the end of the cathode 
cavity. 
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To compare these data with previous measurements of 

enclosed current profiles at the higher currents of 0.9 and 

4.7 kA, the enclosed current data at the higher currents were 

also reduced to yield surface current density distributions. 

Figure 19 shows the resulting current density profiles for a 

—3 

current of 0.9 IcA and mass flows from 5 x 10 to 8 g/sec. 
Here, just as at 0.25 kA, the current density peak moves fur- 
ther upstream from the cavity and and the current attachment 

becomes more diffuse as the mass flow is reduced from 8 to 
—2 

6.6 X 10 g/sec. Again, still further reduction in mass 
flow, from 6.6 X 10 to 5 x 10 g/sec, causes the peak in 
current density to move downstream towards the cavity end. 
However, the maximum penetration at the higher current is only 
0.9 cm, compared to a maximum penetration of l.J cm at a 
current of 0.25 kA. 

Figure 20 shows similar surface current density distribu- 
tion for a current of 4.7 kA, At this current, unlike at the 
lower currents of 0.9 and 0.25 kA, the current density distri- 
bution for different mass flows are observed to be very simi- 
lar. The peak current density occurs somewhere between zero 
and 0.2 cm from the cavity end, and in all cases, the current 
density falls to less than 10% of its peak value within 
0.9 cm. 

Figure 21 stimraarizes the current density measurements by 
graphing the active zone length, defined as the axial distance 
from the cavity end to the peak in current density, against 
the cathode mass flow, with discharge current as a parameter. 
The surface current density distributions are seen to be most 
sensitive to changes in mass flow rate at the lowest current 
of 0.25 kA, becoming least sensitive at the highest current of 

4.7 kA. The maximum penetration of current into the hallow 
cathode cavity decreases as the current is increased, from one 
cathode diameter at 0.25 kA to approximately one- tenth of the 
cathode diameter at 4.7 kA. 
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It is important to note that this decrease of the current 

penetration into the cavity with increasing current proceeds 

despite the quadratically increasing x B force inside the 

cavity, which acts in the upstream direction opposing the gas 
146 

flow. To determine whether this dependence continues to 

even larger currents, the distribution of current was measured 

for a total current of 16.7 kA, where the electromagnetic 

force is known to dominate the acceleration process at these 
138 

mass flows. Figure 22 shows surface current density distri- 

butions for increasing currents of 0.25, 0.9, 4.7 and 16.7 kA, 
each drawn normalized to its peak value. It is observed that 
the trend continues to the higher currents, with the 16.7 kA 
current penetrating even less than the 4,7 kA current, 

A more thorough understanding of the interaction between 
the electromagnetic forces and the gasdynamic forces could be 
obtained by a detailed analysis of the fluid flow processes 
occurring in the hollow cathode. However, such an analysis is 
complicated for two reasons. First, the initial cold flow 
field inside the hollow cathode involves oblicpae, curved shock 
waves whose patterns are difficult to compute. Second, with 
the initiation of the discharge, the cold gas flow field is 
perturbed by electron emission from the cathode wall, electro- 
magnetic body forces, and heat addition via collisional and 
radiative processes. 

Rather than attempt such an analysis, an experiment was 
performed to discriminate between the static and dynamic pro- 
cesses inside the hollow cathode. In this experiment, the 
mass flow to the hollow cathode was completely shut off, and 

the discharge vessel was prefilled to various static pressures 
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of argon from 5 x 10 to 100 torr. At each ambient gas 
pressure, a surface current distribution was obtained for the 
fixed total current of 0.25 kA. These distributions are shown 
in Fig, 23. 
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Comparison of these data with Fig, 18 shows that the dis- 
tributions obtained with various static prefills are closely 
similar to those obtained with mass flow through the hollow 
cathode. As the airibient pressure in the cathode is reduced 
from 100 torr to 5 torr, the surface current density distribu- 
tion becomes more diffuse and the current density peak moves 
upstream to about one cathode diameter from the end of the 
cavity. Further reduction of the pressure from 5 torr to 
5 X 10 torr causes the peak in the current distribution to 
move downstream towards the end of the cavity. 

The effect of discharge current on the surface current 
density distributions at a fixed ambient pressure of 5 torr 
was also investigated. Figure 24 shows the surface current 
density distributions, normalized to the peak current density 
in each case, for currents of 0.25, 0.9, 4,7 and 16.7 1<A, all 
for a fixed ambient pressure of 5 torr. It is observed that 
just as at a mass flow of 0.2 g/sec, Fig, 22, the effect of 
increased current is to move the current distributions further 
downstream towards the end of the cathode. 

The similarity between the distributions of current for 
various mass flows or ambient pressures at fixed current. 

Figs. 18 and 23, and f-r various currents at fixed mass flow 
or ambient pressure. Pigs. 22 and 24, suggests that the dynam- 
ics of the gas flow in the hollow cathode do not play a major 
role in establishing the surface current density patterns. 
Further experiments are planned to aid in discriminating 
between the static and dynamic components of the injected 
mass flow. 

In all of the experiments described above, the variations 
of the surface current density distribution with mass flow, 
current and ambient pressure were determined by traversing 
axially a magnetic probe at a fixed radial position. In order 
to obtain a more comprehensive picture of the energy deposi- 
tion patterns inside the hollow cathode cavity, these surface 
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current density profiles were supplemented by magnetic field 
measurements throughout the cavity and by complete maps of 
plasma floating potential. The distributions were obtained 
for mass flows of 0.4 g/sec and 0,1 g/sec, both at a fixed 
current of 0.25 3cA, Figure 18 shows that these two mass 
flows represent extremes in the penetration of the current 
into the hollow cathode cavity at 0.25 hA. 

In both cases, the plasma floating potentials were meas- 
ured using a single Langmuir probe consisting of an insultated 
0.25-mm-diameter tungsten wire of which only the face is ex- 
posed. The probe output is connected through a 100 M2 P-6013A 
voltage probe to one input of a differential amplifier at the 
oscilloscope. The other input is the cathode potential, meas- 
ured relative to the anode ground with an idential P-6013A 
probe. Thus, the oscilloscope displays the floating potential 
relative to the cathode surface, and since electron temperature 
effects depress the floating potential below the true plasma 
potential, this differential signal is the minimum potential 
drop between the cathode and the probe tip. 

Figure 25a shows, on a cross-sectional view of cathode 
HC-XII, contours of constant current and constant floating po- 
tential for a mass flow of 0,4 g/sec and a current of 0.25 3cA. 
The lower half of the figure. Fig. 25b, shows a similar map 
of current and floating potential contours for the lower mass 
flow of 0.1 g/sec, at the same current. 

For both mass flows, the equipotential lines are approxi- 
mately radial in the bulk of the plasma, i.e. there is a neg- 
ligible radial field in the volume of the cavity. In addition, 
the potentials show a weak axial electric field of less than 
5 V/cm. Since the cathode itself is the zero volt equipoten- 
tial in both cases, all the radial equipotentials must bend 
parallel to the cathode surface soi'-^ /here between the surface 
(R = 0.95 cm} and the largest radial probe position (R = 0.8 cm). 
If the potentials all bend within the Debye sheath separating 
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the surface from the quasi neutral plasma, extremely high 
radial electric fields, of order lO'V/cm, would be produced at 
the surface of the hollow cathode. 

Prom the intersecting grid of current and potential lines, 
the power density, j . ,E, in any volume segment of the plasma 
can be easily determined. The results of this computation are 
graphed in Pig. 26a and b which show power density as a func- 
tion of axial position inside the cavity, with radius as a 
parameter. For both mass flows, it is observed that the power 
density profile at a radius of 0.8 cm is qualitatively similar 
to the surface current density profile measured at this same 
radial position (Pig. 18) . Furthermore, in both cases the 
power density profile retains its form well into the plasma, 
up to approximately half the radius. Thus, the detailed maps 
of current and potential inside the cathode cavity corroborate 
the significant difference in current penetration into the 
cavity inferred from surface current density measurements. 

B . Summary 

Current distributions measured for a range of mass flows 
at a total current of 0.25 kA display trends similar to previ- 
ous measurements at 0.9 and 4.7 3cA. At a fixed current, the 
penetration of the current into the cavity initially increases 
as the mass flow is decreased, but then decreases as mass flow 
is decreased further. The maximum penetration for a given 
mass flow increases as current is decreased from 16.7 to 
0.25 kA, reaching approximately one cathode diameter at 0.25 kA. 
Floating potential maps show that the current penetration is 
accompanied by a weak axial electric field of less than 5 V/cm 
and a radial field which is negligible except near the cavity 
walls. Current distributions for various static prefills with 
no mass flow injected through' the hollow cathode were found to 
replicate those distributions measured with only injected 
flows, yielding the important result that the dynamic head of 
the injected flow does not play an essential role in the 
hollow cathode emission process. 
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APPENDIX A: Semi-annual Statement of Expenditures 


PULSED ELECTRO.MAGNETIC GAS ACCELERATION 
mSA NGL 31-OOl-OOS 

1 January 1975 to June 30, 1975 


DIRECT COSTS 


I. 


II. 

III. 

IV. 

V. 


Salaries 

Professional 
Tecltnicians 
Students 

Supporting Staff 

Employee Benefits ( 22 ^ 5 %) 

Materials and Services 
Travel 
Tuition 

TOTAL DIRECT COSTS 


$ 14,313 
11,765 
3,050 
1,835 


$ 30,963 
6,280 
7,194 
1,477 
2,400 
$ 48,314 


INDIRECT COSTS 


26,783 


VI. Overhead ( 86 ^ 2 %) 


TOTAL 


$ 75,097 



